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北京 地 区 重力 波 活动 及 其 波谱 的 季节 分 布 特性 研究 
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摘 要 利用 激光 雷达 对 北京 地 区 上 空 Na 层 进行 持续 观测 , 通过 连续 三 年 累积 的 夜间 观测 数据 对 北京 地 区 重力 
波 活动 及 其 波谱 进行 研究 . 根据 重力 波 的 线性 理论 计算 , 得 到 北京 地 区 上 空 的 大 气 密度 扰动 规律 、 空 间 功率 谱 和 时 
间 频 率 谱 . 通过 选择 重力 波 波 长 在 1~8km, 具有 特定 波长 以 及 特定 周期 为 60, 45, 25 min 的 重力 波 活动 辅助 研究 
重力 波 的 季节 变化 规律 , 结果 表明 北京 地 区 重力 波 大 气 密 度 扰 动 具有 夏季 大 、 冬 季 小 的 活动 规律 . 结合 波源 与 背景 
风 场 的 季节 性 变化 规律 , 分 析 得 出 北京 上 空 重 力 波 活动 季节 性 变化 的 主要 原因 为 青藏 高 原 地 形 和 对 流 因 素 与 我 国 
北方 地 区 季节 性 背景 风 场 共同 作用 的 结 
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over Beijing are studied with 3 years of sodium lidar observations. From the linear layer density re- 
sponse to gravity wave forcing, the lidar data were analyzed to get the atmospheric density perturba- 
tions and their spectra. The atmospheric density perturbation, density variance for fluctuations with 
vertical scales between 2 and 10km, and amplitudes of density perturbation spectra at m = 27/8, 
27/4, m = 27 /1.5, m = 27/1, and w = 27/60, w = 27/40, 27/25 all exhibit large nightly variability 
as well as large seasonal variations, with the semiannual maxima occurring near the solstice. The 
mean RMS atmospheric density perturbation over Beijing are 5.8%, which are obviously larger in 
summer than that in winter and the maxima occur near the solstice. The m spectra show power 
law shapes, and their range of variation is between —1.97 and —3.67 with an annual mean value of 
—3.02, and w spectra is between —1.06 and —2.08 with an annual mean value of —1.92, respectively. 
It is concluded that the reaction of the Qinghai-Tibet Plateau and the background wind may be the 


main reason of the gravity wave behaviors at Beijing. 


2015, 35(4) 
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0 引言 


作为 中 高 层 大 气 中 最 重要 的 动力 学 过 程 之 一 , 近 
年 来 重力 波 在 中 高 层 大 气 研究 领域 得 到 广泛 关注 . 
对 重力 波 而 言 , 由 于 其 可 将 能 量 和 动量 从 波源 传输 到 
中 高 层 大 气 , 对 中 高 层 大 气 的 循环 、 结 构 、 组 成 都 有 
重大 影响 , 有 助 于 深入 理解 大 气 的 大 尺度 循环 动力 学 
行为 .研究 重力 波 的 产生 及 其 在 耗 散 和 传播 过 程 中 
的 动力 学 规律 , 是 当今 中 高 层 大 气 研究 的 重要 内 容 ， 
对 于 深入 研究 全 球 大 气 能 量 收 支 以 及 能 量 和 动量 平 

激光 雷达 是 研究 重力 波 的 有 效 手 段 , 通过 地 面 发 
射 激光 对 大 气 进行 主动 探测 , 近年 来 在 重力 波 研究 领 
域 得 到 了 广泛 应 用 . She 等 是 利用 测 温 激光 雷达 对 美 
国 Colorado 地 区 中 层 顶 区 的 重力 波 活动 展开 了 观测 
研究 ; Gardner 等 P-5 基于 CEDAR 观测 了 北美 地 区 
重力 波 活 动 , 提取 出 重力 波 活 动 数 据 , 并 总 结 了 Ur- 
bana 地 区 的 重力 波 活 动 行为 ; Yang 等 是 基于 巴西 
激光 雷达 观测 数据 报道 了 南半球 重力 波 季 节 变 化 具 
有 与 北半球 完全 不 同 的 活动 规律 等 . 目前 普遍 认为 ， 
重力 波 的 产生 和 传播 与 观测 地 所 处 地 理 位 置 有 着 极 
其 密切 的 关系 .为 更 清楚 地 研究 重力 波 的 产生 及 其 
传播 特性 , 有 必要 在 全 世界 多 地 范围 内 展开 联合 观测 
研究 . 同时 , 持续 进行 长 时 间 的 观测 数据 积累 有 助 于 
获得 更 加 丰富 的 重力 波 参数 统计 信息 , 从 而 更 好 地 研 
究 其 变化 规律 . 这 对 深入 理解 重力 波 的 行为 特性 是 十 


Lidar, Na layer, Gravity wave, Seasonal variations, Topography of Qinghai-Tibet 


本 文 利用 子午 工程 钠 荧 光 激 光 雷 达 2009—2012 
年 在 中 国 北京 延庆 观测 台 站 获取 的 Na 层 观 测 数 据 ， 
对 北京 地 区 重力 波 活 动 规律 进行 了 研究 ， 通 过 分 析 
北京 地 区 上 空 大 气 密度 扰动 、 空 间 功 率 谱 、 时 间 频 率 
谱 以 及 特定 波长 和 频率 的 重力 波 活 动 行为 , 得 出 北京 
地 区 大 气 密度 扰动 的 季节 变化 规律 , 并 结合 中 国 地 形 
以 及 背景 风 场 的 季节 性 变化 , 分析 了 北京 上 空 重 力 波 


1 数据 与 讨论 


1.1 实验 设备 

子午 工程 激光 雷达 北京 延庆 人 台 站 的 观测 设备 为 
双 波 长 激光 雷达 , 波长 分 别 为 532nm 和 589nm. 所 
用 的 钠 荧光 激光 雷达 数据 采用 染料 激光 名 进行 激发 ， 
产生 589nm 的 激光 探测 光束 (能 量 30mJ), 重复 频 
率 为 30 Hz. 接收 望远镜 参数 为 : 主 镜 直 径 1000mm, 
发 射 激光 光束 距离 接收 望远镜 中 心 6m, 望远镜 视 场 
通常 设置 为 2 mrad, 接收 信号 通过 带宽 1nm 的 带 通 
光学 滤 光 片 、 致 冷光 电 倍 增 管 、 前 置 放 大 器 和 时 间 分 
辨 的 光子 计数 融 进 行 信号 探测 与 放大 , 最 终 由 计算 机 

每 个 激光 雷达 文件 由 5000 个 激光 脉冲 累积 , 高 
度 分 辩 率 为 96 m, 时 间 分 辨 率 约 为 167s. Na 层 密 度 
可 由 这 些 原始 激光 雷达 回 波 数据 得 到 M. 
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1.2 ”数据 处 理 

利用 钠 激 光 雷 达 进 行 Na 层 数 据 提取 是 通过 Na 
层 激 光 雷 达 研 究 重 力 波 的 关键 .本 文采 用 Yang 
等 5s 9] 改进 的 Na 层 数 据 提取 方法 ,通过 对 当天 Na 
密度 进行 处 理 ,求解 得 到 本 底 Na 层 密度 ,进而 计 
算 Na 层 密 度 扰 动 ， Na 层 扰 动 与 大 气 密度 扰动 关 
系 为 

rs(z,t) = -一 | ae 

A, rs(z,t) 为 Na 层 密度 扰动 , ra(z,t) 为 大 气 密度 
TB), > 和 上 分别 代 表 实 时 观测 到 的 Na 层 高 度 及 观 
测 时 间 , co 为 未 扰动 Na 层 的 平均 厚度 , 约 为 4.5km， 
zo 为 钠 层 中 心 高 度 . 通常 在 距离 地 表 94 士 2km 高 度 
Mh, Na 层 密度 较 大 且 随 重力 波 的 传播 变化 剧烈 . 因 
此 , 在 提取 重力 波 参 数 时 将 这 一 部 分 Na 层 密度 数据 
EF. 最 终 大 气 密 度 扰 动 可 利用 

Ay id 
TS 9H fa) 
计算 . AF, f(z) = —po(2)/po(2), 即 本 底 Na R 
度 与 本 底 钠 层 密 度 微分 之 比 . 

对 大 气 密度 扰动 进行 全 天 时 和 全 空间 平均 化 处 
H, 得 到 大 气 密度 的 均 方 根 值 , 即 


ralz t): 


(1) 


ralz, t) = rs(z,t) 


(2) 


P 1 | 三 
T(t) = 
( 3 )) TL! zo — L/2 
(3) 


式 中 , L! 为 舍弃 位 于 942km 处 非 高 斯 分 布 Na 层 
厚度 后 的 总 宽度 值 , 7 为 当天 的 总 观测 时 间 . 通过 子 


r2(z,t)dtdz. 
t—T/2 
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午 工程 北京 延庆 站 累积 的 2009—2012 年 共 173 RA 
效 观测 数据 , 经 上 述 过 程 反 演 并 平均 化 求 得 北京 地 区 
大 气 密 度 扰 动 季节 分 布 值 ， 将 所 有 统计 数据 进行 最 
小 均 方 误差 (MMSE) 拟 合 后 得 出 北京 地 区 大 气 密度 
扰动 的 年 度 季 节 变 化 规律 , 如 图 1 所 示 . 

图 1 给 出 的 是 北京 地 区 大 气 密度 扰动 年 度 均 方 
根 分 布 , 大 气 扰 动 拟 合 参 数列 于 表 1. 实验 数据 表 
H, 年 度 大 气 密度 扰动 平均 信 为 5.83%, 并 且 夏 季 的 
大 气 密度 扰动 显著 高 于 其 他 季节 . 全 年 扰动 贡献 值 
为 0.98%, 优 于 半年 扰动 贡献 0.72%, 分 别 占 全 年 平 


(107) 


RMS atmosphere density pertubation 


120 


180 240 
day number 
图 1 北京 地 区 大 气 密度 扰动 年 度 均 方 根 分 布 ( 圆 点 代表 

气 密度 扰动 的 均 方 根 数据 , 黑色 曲线 表示 经 
MMSE 拟 合 的 大 气 扰动 谱 线 ) 


Fig.1 Seasonal variations of the RMS atmospheric 


density perturbations (Circles represent the RMS 
atmospheric density perturbations and the solid curve 
represents the MMSE fitting of the RMS 


atmospheric density perturbations) 


表 1 重力 波 拟 合 参数 年 度 及 半年 度 贡 献 值 


Table 1 Annual and semiannual components of gravity wave parameters 


(72) /2/(%) 


Ao Ai A2 

5.84 0.98 0.72 

Em/(%)? 5.13 1.85 1.37 

Fa (ms) /(m-cycle~*) 4.47 1.69 1.21 
Fa (ma) /(m-cycle~*) 1.20 0.58 0.34 
F,(m1.5)/(m-cycle~*) 0.61 0.376 0.207 
Fa (m1) /(m-cycle~*) 0.013 0.008 0.0045 
F, (woo) /(s-cycle7') 1.34 0.786 0.67 
F, (was)/(s-cycle') 0.761 0.51 0.37 
F, (wis)/(s-cycle~*) 0.124 0.093 0.052 


dı də AyfAy Az/Ao d-d 
181 188 0.16 0.12 ae 
198 179 0.36 0.26 19 
193 187 0.39 0.27 6 
197 187 0.48 0.28 10 
193 188 0.61 0.34 5 
190 194 0.60 0.33 一 4 
174 193 0.59 0.50 —19 
180 188 0.67 0.49 -8 
179 186 0.75 0.41 z 


注 重力 波 参 数 方程 y = Ao + Aı cos[27/365(d — d1)] + A2 cos[4r /365(d — d2)]. 
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均 扰 动 的 16% 和 12%. 大 气 密度 扰动 的 最 低 值 出 现 
在 春分 和 秋分 日 前 后 . 总 体 呈 现 夏 季 大 , 冬季 小 的 分 

与 文献 对 比 可 知 , 本 文 结果 与 Senft 等 报道 的 美 
E| Urbana 地 区 (56°W, 45°N) 大 气 密度 扰动 的 季节 
规律 类 似 站 .其 中 Urbana 地 区 大 气 密度 扰动 的 平 
均值 为 5.6%, 与 本 文 数值 十 分 接近 并且 都 是 在 夏 
季 扰 动 幅 度 最 大 , 而 春秋 分 时 出 现 扰动 极 小 值 . 


2 重力 波 参数 拟 合 值 


研究 重力 波 场 能 量 在 各 个 垂直 波 数 或 频率 上 的 
相对 分 布 , 可 以 计算 出 重力 波 的 垂直 波 数 谱 和 时 间 频 
率 谱 . 其 中 , 垂直 波 数 谱 (功率 谱 ) 可 采用 自 相关 函数 
法 进行 计算 , 有 

Falm) -| B,(s, 0)e"Sds = 


| (ralz, t)ra(z — s,t))e'™sds. 


AF, 大 气 密度 扰动 的 自 相关 函数 等 于 (rs(z,t)ra(z 一 
s,t)), 可 以 用 Ba(s,0) 表示 , 其 中 为 两 个 采样 的 
高 度 差 值 ， 根据 Tsuda 等 DO 的 处 理 方法 ， 先 对 
大 气 密度 扰动 进行 预 白 噪声 化 处 理 ， 预 白 噪声 化 函 
数 y(z,t) = ra(z,t) — 0.95ra(z — Az, t); 再 对 其 进行 
加 Hamming 窗 修正 后 , 经 傅 里 叶 变换 即 可 得 到 北京 
地 区 每 一 单 天 事件 的 波 数 谱 及 (mm)， 同 时 计算 出 重 
力 波 在 各 个 垂直 波长 分 布 的 谱 斜 率 p 及 其 年 均 统计 
分 布 ， 通 过 对 所 有 单 天 事件 的 波 数 谱 曲线 进行 平均 
化 处 理 , 得 到 年 平均 的 谱 斜 率 p = 一 3.02, 如 图 2(a) 
所 示 . 这 与 Dewan 和 Good! 提出 的 线性 不 稳定 理 
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论 和 Hines? 提出 的 多 普 勒 扩散 理论 预测 的 p = -3 
符合 得 非常 好 . 垂直 波 数 谱 斜率 的 季节 性 统计 分 布 规 
律 如 图 2(b) 所 示 . 其 中 最 缓 的 斜率 为 -1.96, 最 陡 
的 斜率 为 -3.67. Senft 和 Gardner?) 提出 的 斜率 值 
在 —2.20~—3.55 内 , 变化 范围 比 本 文 小 . Senft 等 [9 
利用 一 个 月 的 数据 研究 了 最 组 斜率 为 -2.1 的 波 事 
件 , 其 结果 相对 接近 本 文 , 而 其 最 大 值 -3.5 低 于 北 
京 地 区 的 观测 值 . 其 他 例如 Collins 等 H4 得 到 的 南 
极地 区 功率 谱 和 斜率 范围 为 一 2.0~ 一 3.5, 也 与 本 文 结 果 
比较 接近 . 

大 气 密 度 方差 Em 是 反映 大 气 重力 波 扰动 的 一 
个 有 效 参 量 , 波长 低 于 10km 的 重力 波 可 以 反映 出 宽 
为 20km Na 层 的 可 靠 信 息 . 本 文 定义 的 Em 是 通过 
选取 波长 范围 在 1~10km, 并 对 Fam) 进行 积分 得 
到 的 , 如 图 3 所 示 , 北京 地 区 Em 最 大 值 为 24x10-4 
最 小 值 0.11x10-4, 全 年 平均 值 为 4.7x10-4 最 
大 值 发 生 在 夏季 的 夏至 点 时 ， 半 年 度 和 年 度 贡 献 值 
分 别 为 1.85x10-? 和 1.37x10-?， 其 全 年 及 半年 变 
化 值 约 占 全 年 波动 的 36% 和 26%, 均 小 于 Senft 
和 Gardner?! 报告 的 45% 和 44%. 如 图 3 Bras, db 
京 地 区 大 气 密 度 扰动 的 Es 季节 分 布 规律 与 均 方 根 
值 分 布 规律 大 体 一 致 , 夏季 的 大 气 密度 扰动 显著 高 于 
其 他 季节 , 最 低 值 出 现在 春分 及 秋分 日 前 后 . 其 中 秋 
分 日 时 的 扰动 更 强 一 些 . 

通过 对 波长 为 1, 1.5, 4, 8km 的 重力 波 活动 强度 
进行 研究 , 给 出 了 其 年 度 分 布 规律 ( 见 图 4). 上 述 波 
长 范围 在 1~8km 的 特征 波长 垂直 波 数 谱 的 季节 变 
化 特性 与 大 气 密 度 扰动 的 季节 变化 特性 一 致 . 如 图 4 
所 示 . 分 别 用 ms = 27/8, ma = 27/4, mi 5 = 27/1.5, 
mi 二 27/1 来 代表 . 由 图 4 和 表 1 可 以 看 到 , ms 的 年 


mean: -3.00 min:-1.95 max: -3.67 


-2.0 上 (b) 
| 
Be “°F a The 
sa ; fae ja 
S2 -3.0 Lee” ey Supt i Tel Dems A 
SE | ee yt 
_ 2 ke . 
ES -35i 
5 ` . 
2 L 

-4.0 L aT 

0 60 120 180 240 300 360 


day number 


图 2 垂直 波 数 谱 的 年 度 平均 谱 (a) 和 垂直 波 数 谱 斜 率 (b) 的 统计 分 布 


Fig.2 Distribution of annual mean of vertical wave number spectra (a) and the spectrum slopes (b) 
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度 平 均值 为 4.47mcycle-1, 远 小 于 Senft 等 外 提 
出 的 12.5 mcycle™, 仅 为 其 36%. 年 度 和 半年 度 振 
幅 分 别 占 全 年 平均 的 39%, 27%. ma 的 年 度 平 均值 
是 1.20m'cycle ,全 年 变化 的 相对 振幅 是 48%, 半年 
值 为 28%. mis, m 为 小 波长 重力 波 扰动 , 全 年 平均 
值 分 别 为 0.61 和 0.013m'cycle ,其 全 年 相对 振幅 
很 高 , 均 达 到 60% 以 上 , 远 高 于 其 余 尺度 波 长 扰动 比 
例 . 值得 注意 的 是 , 所 有 波长 的 波动 行为 均 是 全 年 振 
幅 比 例 高 于 半年 振幅 比 , 其 相位 也 出 现在 一 年 的 中 间 
时 , 即 夏至 点 附近 . 这 一 结论 与 前 文 所 述 大 气 密度 扰 
动 、 波 数 谱 扰 动 方差 的 季节 分 布 规律 均 符合 得 很 好 . 
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时 间 频 率 谱 
Falw) = ie (ra(z,W)ralz — s, w))e* dz 


也 是 研究 重力 波 扰动 的 重要 参量 ， 与 大 气 扰动 的 垂 
直 功 率 谱 相 似 , 在 对 大 气 密 度 扰动 进行 漂白 及 加 徐 
处 理 后 , 通过 侍 里 叶 变 换 得 到 所 有 单 天 事件 的 频率 
谱 Faw). 通过 计算 谱 和 斜率 4 及 其 年 均 统计 分 布 发 
M, g 的 分 布 范 围 为 一 1.06~ 一 2.80, 平均 值 为 一 1.92 
( 见 图 5). Gardner 等 P) 得 出 的 北美 地 区 当地 值 
为 —1.92, Beatty 等 15 观测 到 的 Arecibo 地 区 值 
为 -1.85, 本 文 结果 与 上 述 结果 均 较 为 接近 , 并 接近 
于 理论 预测 值 —5/3. 

通过 对 波动 周期 分 别 为 60, 40, 25min 的 重力 
波 活动 的 季节 变化 研究 , 给 出 了 其 年 度 分 布 规律 如 
图 6 Bras. ”由 图 6 AIR 1 可 以 看 到 ，weo，wao， 
wos 的 年 度 平均 值 分 别 为 1.34, 0.76, 0.12s-cycle7t, 
其 波动 行为 均 是 全 年 振幅 比例 要 高 于 半年 振幅 比 ， 
相位 也 出 现在 夏至 点 附近 . ”这 说 明 上 述 频 率 范围 
在 27/60~27/25 的 特定 波长 的 垂直 波 数 谱 季节 变化 
特性 与 大 气 密度 扰动 的 季节 变化 特性 一 致 . 

研究 结果 显示 ， 北 京 地 区 的 重力 波 季节 变化 
特性 ， 例 如 大 气 密度 扰动 、 密 度 扰 动 方 差 、8km 
和 60min 以 下 范围 的 垂直 功率 谱 及 时 间 频 率 谱 ， 
均 呈 现 夏季 大 而 冬季 小 的 季节 变化 规律 . 这 一 结果 
与 Gong 等 9 报道 的 重力 波 事件 的 出 现 频率 相 吻 
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图 4 不 同 波长 重力 波 垂直 波 数 的 季节 变化 规律 


Fig.4 Seasonal vibrations of vertical wave number spectral amplitudes 
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图 5 时 间 频 率 谱 的 年 度 平均 谱 (a) 和 时 间 频 率 谱 和 斜率 (b) 统计 分 布 
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Distribution of annual mean of temporal frequency spectra (a) and the spectrum slopes (b) 
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图 6 不 同 频率 的 重力 波 垂直 频率 季节 变化 分 布 


Fig.6 Seasonal vibrations of temporal frequency spectral amplitudes 


合 . 文献 研究 表明 , 地 形 和 对 流 作 用 是 产生 重力 波 的 
主要 因素 07-25], Wan 等 5 和 Xu 等 PO) 在 基于 广 
泛 的 HF Doppler 阵列 观测 基础 上 报道 了 青藏 高 原 上 
的 对 流 和 地 形 波 是 中 国 中 部 地 区 上 空 大 气 中 尺度 重 
力 波 的 主要 波源 , 而 中 层 大 气 中 观测 到 的 重力 波 主要 
HEA 20km 以 下 的 低层 大 气 , 并 且 大 气 对 流 是 重力 波 
中 常见 的 主要 波源 . 北京 地 区 地 处 中 国 北方 , 对 流 层 
温度 夏季 较 高 而 冬季 较 低 . 因此 在 夏季 , 对 流 强度 和 
频 度 都 是 一 年 中 最 大 的 ; 而 在 冬季 , 对 流 强 度 和 频 度 
则 是 一 年 中 最 小 的 . 也 就 是 说 , 波源 的 强度 和 频 度 在 
一 年 中 夏季 最 大 而 冬季 最 小 . 

另外 , 根据 Gongl'9 基于 HWM93* 描绘 出 的 
中 国 北方 地 区 水 平 风 场 模型 OO) 并 辅助 结合 Zhang 


*http://modelweb.gsfc.nasa.gov/ 


和 YiR9 基于 Radiosonde 描绘 的 北京 地 区 对 流 层 风 
场 的 变化 规律 , 可 以 看 出 中 国 北方 背景 风 场 的 季节 变 
化 , 尤其 是 水 平 纬 向 风 场 的 季节 变化 , 具有 明显 的 半 
年 对 称 性 (水 平 经 向 风 比 水 平 纬 向 风 小 得 多 , 水 平 经 
向 风 对 重力 波 活动 的 影响 比较 弱 ， 因 此 在 分 析 背 景 
风 场 的 影响 时 , 只 考虑 了 纬 向 风 的 作用 ). 在 春季 和 夏 
季 , 中 国 北方 上 空 纬 向 风 为 西向 风 , 这 一 季节 的 背景 
大 气 对 向 东北 方向 传播 的 重力 波 的 吸收 、 耗 散 作 用 
较 弱 ; 而 秋季 和 冬季 纬 向 风 为 东 向 风 , 由 于 Doppler 
效应 , 此 时 背景 大 气 对 重力 波 的 滤波 效应 最 强 . 由 此 
可 知 , 在 背景 风 场 方 向 与 重力 波 水 平 传播 方向 一 致 的 
季节 , 由 于 背景 大 气 的 吸收 和 耗 散 效应 使 得 大 部 分 重 
力 波 消 失 于 大 气 中 , 雷达 观测 到 的 大 气 重力 波 扰动 较 
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小 ; 在 背景 风 场 方向 与 重力 波 传播 方向 相反 的 季 市 ， 
背景 大 气 对 重力 波 的 吸收 与 耗 散 作用 较 小 , 重力 波 引 
起 的 大 气 密度 扰动 很 大 ， 因 此 重力 波 扰动 在 夏季 出 
现 最 大 值 , 而 冬季 观测 到 的 重力 波 大 气 扰动 则 相对 较 
弱 . 从 冬季 到 夏季 , 由 于 对 流 引 起 的 波源 强度 和 频 度 
随 着 时 间 变 化 而 逐渐 增强 , 地 面 附近 的 大 气温 度 也 随 
时 间 变 化 而 逐渐 升 高 ,导致 重力 波 的 扰动 强度 增 大 . 
同时 , 中 国 上 空 纬 向 风 逐 渐 由 东 向 风 变 为 西向 风 , 背 
景 大 气 对 传播 的 重力 波 的 阻碍 作用 逐渐 减弱 ; 相反 ， 
从 夏季 到 冬季 , 由 于 对 流 引 起 的 波源 强度 和 频 度 随 时 
间 变 化 而 逐渐 减弱 , 地 面 附近 的 大 气温 度 也 随 着 时 间 
变化 而 逐渐 降低 , 大 气 扰动 的 强度 也 因 之 减弱 . 与 此 
同时 , 中 国 上 空 的 纬 向 风 逐 渐 由 西向 风 变 为 东 向 风 ， 
背景 大 气 对 传播 的 重力 波 的 阻碍 作用 逐渐 增强 , 导致 
秋冬 季节 观测 到 的 大 气 密度 扰动 相对 较 弱 . 

基于 以 上 分 析 认 为 , 北京 地 区 上 空 重 力 波 活动 季 
节 性 变化 的 主要 原因 是 青藏 高 原 的 地 形 及 对 流 因素 
与 背景 风 场 的 季节 性 变化 共同 作用 的 结果 从 冬季 
到 夏季 , 雷达 观测 的 重力 波 扰动 强度 不 断 增 大 ; 从 夏 
季 到 冬季 , 可 观测 到 的 重力 波 活动 强度 越 来 越 弱 . 也 
就 是 说 , 从 冬季 到 夏季 重力 波 活动 性 逐渐 增强 , 从 夏 
季 到 冬季 , 重力 波 活 动 性 逐渐 减弱 , 具有 明显 的 夏季 
极 大 冬季 小 的 季节 变化 性 . 另外 , 本 文 的 重力 波 研 究 
数据 受 激光 雷达 自身 特点 限制 , 所 包含 的 信息 不 够 全 
H, 探索 重力 波 季节 性 变化 的 根本 原因 还 需要 更 多 观 
测 手段 , 这 有 待 进一步 深入 研究 . 


利用 激光 雷达 技术 对 北京 地 区 上 空 Na 层 进行 
观测 , 积累 连续 三 年 的 夜间 观测 数据 对 重力 波 活动 
规律 进行 研究 ， 根 据 重力 波 的 线性 理论 ， 计算 得 到 
了 北京 地 区 上 空 大 气 密度 扰动 、 空 间 功 率 谱 及 时 间 
频率 谱 . 通过 选择 波长 为 8 4 1.5, 1 km 以 及 波动 
周期 为 60, 40, 25min 的 重力 波 辅助 研究 大 气 密度 
扰动 的 季节 变化 规律 ,得 到 大 气 密 度 扰动 具有 夏季 
K, 冬季 小 , 峰值 约 出 现在 夏至 时 的 结论 , 平均 大 气 
密度 扰动 的 均 方 根 值 为 5.8%. 通过 加 Hamming 窗 
滤波 平滑 处 理 , 拟 合 得 到 垂直 功率 谱 p 值 变化 范围 
为 一 1.97~ 一 3.67, 年 平均 值 为 -3.02. 频率 谱 的 谱 斜 
率 q 值 分 布 范围 为 -1.06~ 一 2.80, 平均 值 为 —1.92. 
该 结论 表明 ， 北京 地 区 重力 波 活 动 规律 与 扩散 滤波 
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理论 相符 .根据 中 国 北方 地 区 背景 风 场 季节 变化 与 
水 平 纬 向 风 场 季节 变化 规律 , 认为 青藏 高 原 的 地 形 和 
对 流 因 素 以 及 背景 风 场 的 季节 性 规律 是 导致 北京 地 
区 上 空 重 力 波 活动 季节 性 变化 的 主要 原因 . 


致谢 分 析 数 据 得 到 了 国家 子午 工程 和 北京 延庆 激光 雷 
达 台 站 的 支持 . 
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